Superconducting proximity effect to the block antiferromagnetism in Kj/Fe2 xSe2 
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Recent discovery of superconducting (SC) ternary iron selenides has block antiferromagentic 
(AFM) long range order. Many experiments show possible mesoscopic phase separation of the 
superconductivity and antiferromagnetism, while the neutron experiment reveals a sizable suppres- 
sion of magnetic moment due to the superconductivity indicating a possible phase coexistence. 
Here we propose that the observed suppression of the magnetic moment may be explained due to 
the proximity effect within a phase separation scenario. We use a two-orbital model to study the 
proximity effect on a layer of block AFM state induced by neighboring SC layers via an interlayer 
tunneling mechanism. We argue that the proximity effect in ternary Fe-selenides should be large 
because of the large interlayer coupling and weak electron correlation. The result of our mean 
field theory is compared with the neutron experiments semi-quantitatively. The suppression of the 
magnetic moment due to the SC proximity effect is found to be more pronounced in the d-wave 
superconductivity and may be enhanced by the frustrated structure of the block AFM state. 

PACS numbers: 74.20.Mn, 74.25.Ha, 74.62.En, 74.25.nj 



I. INTRODUCTION 



The recent discovery of high- Tc superconductivity 
in the ternary iron selenides Aj^Fe2_2:Se2 (A=K; Rb; 
Cs;...)i"— has triggered a new surge of interest in study of 
iron-based superconductors (Fe-SC). The fascinating as- 
pect of these material lies in the tunable Fe- vacancies in 
these materials, which substantially modifies the normal- 
state metallic behavior and enhances the transition tem- 
perature Tc to above 30K from 9K for the binary sys- 
tem FeSe at ambient pressure jiisi Particular attention 
has been focused on the vacancy ordered 245 system, 
K0.8Fe1.6Se2, as it introduces a novel magnetic struc- 
ture into the already rich magnetism of Fe-SC. Unlike the 
collinear— or bi-coUinear— AFM order observed in the 
parent compounds of other Fe-SC, the neutron diffraction 
experiment has clearly shown that these materials have 
a block AFM (BAFM) order Meanwhile, the AFM or- 
der with an unprecedentedly large magnetic moment of 
3.31/iB/Fe below the Neel temperature is the largest one 
among all the known parent compounds of Fe-SCfiiii^ 
Moreover, the carrier concentration is extremely low, in- 
dicating the parent compound to be a magnetic insu- 
lator/semiconductor ^i^ii^ in comparison with a metallic 
spin-density-wave (SDW) state of the parent compound 
in other Fe-SC^'^^ 

The relation between the novel magnetism and super- 
conductivity in ternary Fe selenides is currently an in- 
teresting issue under debate. The question is whether 
the superconductivity and the BAFM order are phase 
separated or co-exist in certain region of the phase di- 
agram. The neutron experiment shows the suppression 
of the AFM ordering below SC transition point sug- 
gesting the coexistence. Some other experiments, such 



as two-magnon Raman-scattering, and muon-spin ro- 
tation and relaxation are consistent with this picture. 
On the other hand, the ARPES,i§ NMR M and TEM 
— experiments indicate a mesoscopic phase separation 
between the superconductivity and the insulating AFM 
state. Most recently, Li et al. showed the superconduc- 
tivity and the BAFM orders to occur at different layers 
of the Fe-selenide planes in the STM measurement.^^ 

The vacancy in Fe-selenides is an interesting but com- 
plicated issue. The vacancy in the Fe-selenide carries a 
negative charge since the Fe-ion has a valence of 2+. In 
the equilibrium, we expect the vacancies to repel each 
other at short distance for the Coulomb interaction and 
to attract to each other at a long distance for the elastic 
strain. Such a scenario would be in favor of the phase sep- 
aration to form a vacancy rich and vacancy poor regions 
in the compound. The challenge is then to explain the ob- 
served suppression of magnetic moment of the BAFM due 
to the superconductivity. At the phenomenological level, 
the suppression of magnetism due to superconductivity 
has been reported previously)^ and such phenomenon 
may be explained by Ginzburg-Landau theory.— 

In this paper, we propose that the proximity effect 
of superconductivity to the BAFM in a mesoscopically 
phase separated Fe-selenides may be large to account 
the suppression of the AFM moments observed in neu- 
tron experiment. More specifically, we use a microscopic 
model to study the proximity effect on a layer of the 
BAFM state induced by adjacent SC layer. The prox- 
imity effect in Fe-selenides is expected to be important 
for the two reasons. One is the weaker correlation ef- 
fect, and the other is the larger interlayer hopping am- 
plitude, compared with those in cuprates. Both of them 
may enhance proximity effect on the magnetism from the 



neighboring SC layer. Our model calculations show the 
proximity effect in a mesoscopically phase separated state 
of Fe-selenides may explain various seemingly conflicted 
experiments. 



II. MODEL AND MEAN FIELD THEORY 



AyFe2-a;Se2 is a layered material with FeSe layers sepa- 
rated by alkali atoms, similar to the 122 material in iron 
pnictides family. To investigate the proximity effect to 
the BAFM layer, we consider a single BAFM layer next 
to a SC layer as shown schematically in Fig. [TJ The elec- 
tronic Hamiltonian describing the BAFM layer is given 

by 



BAFM 




H — Hq + Hi, 



(1) 



where Hq describes the electron motion and spin cou- 
plings in the BAFM layer and Hmter describes the cou- 
pling to the neighboring SC layer. We consider a two- 
orbital model to describe iJo, 



<ij>,al3 <<ij>>,al3 
+ Jl Si^a ■ Sj^i3 + J2 Si^Q • Sj^^ 



<ij>' ,af3 



«ij>>' mI3 



where Ci^aa annihilates an electron at site i with orbital 
Oi {dxz and dyz) and spin a, /i is the chemical potential. 
iij,a/3 are the hopping integrals, and < ij > (< ij >') and 
<< ij >> (<< ij >>') denote the intra-block (inter- 
block) nearest (NN) and next nearest neighbor (NNN) 
bonds, respectively [see the upper layer in Fig.[T]. Ji (J2) 
are the exchange coupling constants for NN (NNN) spins 
in the same block, and J{ {J'2) are for the two NN (NNN) 
spins in different blocks. The two-orbital model is a crude 
approximation for electronic structure. However, it may 
be a minimal model to capture some of basic physics 
in examining the proximity effect. The band structure 
around the obtained Fermi energy from the two-orbital 
model is shown in Fig. [21 which is very similar to the 
result obtained in density functional theory. The main 
shortcoming in using the two-orbital model is that the 
magnetic moment is 2/iB at largest, smaller than the ex- 
perimentally measured 3.31/xB. We consider this to be a 
quantitative issue, and will not qualitatively change our 
results. 

We now consider Hinter, the coupling between SC layer 
and BAFM layer. Because of the semiconducting gap of 
the BAFM layer and the SC gap in the SC layer, the 
leading order of the interlayer coupling are the pairing 
hopping between the SC and BAFM layers. According 
to the crystal structures, such a coupling term can be 



FIG. 1: (color online) Schematic diagram of the system in 
our model to study proximity effect to the block AFM state 
(upper layer) induced by superconductivity at the lower layer 
via a pair tunneling process Hinter in Eq. (3). 



expressed as 



Hi, 



/2 

+H.C.), 



(3) 



where ujc is a characteristic energy, and is the interlayer 
hopping integral, the superscript B represents the SC 
layer. With the mean field approximation A^,- 



CL.^Cfa^',<f >' have 



= < 



H^^ter^ E iVr,'jClp^.Clp,^^+H.C.), (4) 
where V;,^ = ^ Y.aa' ^ij,c.a',a5- 

The X 75 vacancy order and the BAFM order lead 
to an enlarged unit cell with eight sites per unit cell. We 
use a mean field theory for the Ising spins in Eq. (2) 
and obtain the Bogoliubov-de Gennes equations in the 
enlarged unit cell 



2^ I H* 



H, 



X exp[zk . (r, - r,)] ( ^^'^'^ ) = f 5'""'" ) '^^^ 

where, the summation of k are over the reduced Brillouin 
zone, and 



(6) 
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Here, JT,intra = Ji (^r.mter = J[) if T = ±x , ±y and 
Jr^zntra = J2 {Jr^znter = J2) if T ±x±y. X and jj denote 
the unit vectors along x and y directions, respectively. 
< S'i+r,/3 > is defined as (n^+^^^^t " 'T'i+r,/3,;)/2. Un,j,a,<7 
Kj,/3,ff)' "«j,a,a (%li,/3,3) are the Bogoliubov quasi- 
particle amplitudes on the j-th site with correspond- 
ing eigenvalues E^. The self-consistent equations of the 
mean fields are 

The magnitude of the magnetic order on the i-th site and 
the induced SC pairing correlation in the BAFM layer are 
defined as, 

k,n 

+<^a,*"!;,,,/5..e^*'-^"^^^"^'^)tanh(^)(8) 

In the calculations, we choose the hopping integrals as 
foUows:^^ Along the y direction, the dxz — dxz NN hop- 
ping integral ti — 0.4 eV and the dyz — dyz NN hopping 
integral t2 = 0.13 eV; they are exchanged to the x di- 
rection; the NN interorbital hoppings are zero; the NNN 
intraorbital hopping integral = —0.25 eV for both d^z 
and dyz orbitals, and the NNN interorbital hopping is 
ti = 0.07 eV. The hopping integral ti is taken as the en- 
ergy unit. We keep Ji : J( : J2 : ^2 = "4 : -1 : 1 : 2^i2i 
The doping level is given by J = n — 2.0. 

III. RESULTS 

To begin with, we present the energy band structure 
at half filling with n = 2 in Fig. [D^a), where Ji = 2.0 is 
so chosen to get a band gap '-^ 500meV being in agree- 
ment with the first principle calculationsi ^^'^^ For the 
electron doping with n = 2.1, the Fermi level crosses an 
energy band around the center of the Brillouin zone [F 
point in Fig.[2jb)], while it intersects with an energy band 
around the zone corner at the hole doping with n = 1.9 
\M point in Fig. djc)]. Although a simple two-orbital 
model is adopted here, both the electron and hole dop- 
ing cases with (5 = 0.1 are qualitatively consistent with 
the first principle calculations.^- In the presence of the 
ordered vacancies and BAFM order, the original two- 
band structures are splitting to sixteen subbands as a 
result of the enlarged unit cell with 8 sites. At half fill- 
ing, 8 lower bands are occupied, i.e., 1/4 electron per one 
subband, while another 8 bands above the Fermi energy 
are unoccupied, resulting in a band gap in Fig. [5] For 



the electron and hole doping with 6 = 0.1, the chemical 
potential crosses one subband which produce the char- 
acteristic features of the Fermi surface and the metallic 
BAFM state. 




r X M r r X M r r x M r 

FIG. 2: (color online) Electronic band structures of Ho given 
by Eq. (2). The parameters are given at the end of section II 
of the text, and Ji = 2.0. (a): at half filling or n = 2.0; (b): 
at electron doping n = 2.1; and (c): at hole doping n = 1.9. 
The color scale indicates the relative spectra weight. 



Motivated by the agreement of the self-consistent 
mean-field solutions with the mentioned first principle 
calculations, we consider now the proximity effect in 
BAFM layer induced by the SC in SC layer. For ex- 
plicit reason, we choose two possible singlet pairing sym- 
metries in the SC layer, i.e., the NNN s±-wave and the 
NN d-wave symmetries with their respective gap func- 
tions As^ — Aq cos{kx) cos{ky) and = Ao[cos(fc^) — 
cos(fcj,)], where the former results in the NNN bond and 
the latter the NN bond couplings in the BAFM layer. 
The interlayer hopping constant tr is assumed to be site 
independent. Fig. |3] displays the moment of the BAFM 
order as a function of the effective tunneling strength 
Vr.ij- At the half filling, both symmetries of the SC or- 
der in the SC layer introduce the decrease of the BAFM 
order and simultaneously induce the SC correlation with 
the same symmetries in the BAFM layer as the tunnel- 
ing strength increases. A main difference between the 
s±- and the d-wave symmetries is the more pronounced 
proximity effect in reducing the moment of the BAFM 
order produced by the c?-wave symmetry as the tunnel- 
ing strength increase, as shown in Figs, ^a.) and[31Jb). 
In the case of electron doping with n = 2.1, where the 
metallic BAFM state results, although the proximity ef- 
fect is more pronounced, the magnetic and the induced 
SC correlation remain the qualitatively unchanged, due 
possibly to the very low total carrier concentration. 

The unique feature of the the effective tunneling in 
the second order is it's temperature dependence via the 
SC pairing Aij^aa' ^a-,cr' , which differs from that in one 
particle tunneling processi^^— The temperature depen- 
dence of the SC pairing parameter is modeled by a 
phenomenological form with A = Aq^I — T/Tc- We 
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FIG. 3: (color online) Block AFM moment and the SC pair- 
ing correlation as functions of the effective tunneling strength 
VT,ij- Black curves are for next nearest neighbor s±-wave 
pairing, and red for nearest neighbor d-wave pairing. Up- 
per panel (a) and (b): n = 2.0 and lower panel (c) and (d): 
n = 2.1. 



present the temperature dependence of the magnetic mo- 
ment in Fig. mja) for the typical choice of the coupling 
constants gr = 2Vr,ijAo = 0.25. As temperature de- 
crease, the magnetic order increases when temperature is 
above Tc, while it decreases when temperature is below 
Tc, resulting in a broad peak around Tc- We note that the 
temperature dependence of the AFM moment is reminis- 
cent of the neutron diffraction and the two-magnon ex- 
periments [Fig. Sl^b)] jiiii^ There is another scenario that 
the competition between the AFM and the SC orders in 
the microscopic coexistence of them may also produce 
the decrease of the AFM moment below Tc- The study 
of such possibility is currently under way and the results 
will be published elsewhere. It is worthwhile to notice 
that the sizable proximity effect relies on the substantial 
interlayer hoping constant t^. Based on the first princi- 
ple calculation, the interlayer hopping tr was estimated 
to have a comparable magnitude with ti possibly due to 
the high values of electron mobility from the intercalated 
alkaline atoms and leads to the highly three dimen- 
sional Fermi surface i^i^^ 

IV. SUMMARY AND DISCUSSIONS 

In summary, we have proposed that various seemingly 
conflict experiments on the phase separation or coexis- 
tence of superconductivity and BAFM may be explained 
within a phase separation scenario by taking into ac- 
count of the proximity effect of superconductivity to the 
neighboring layer of BAFM. We have theoretically stud- 
ied the proximity effect to a BAFM layer induced by 
adjacent SC layers in a simplified two-orbital model for 
Fe-selenides. The proximity effect in reducing the mo- 
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FIG. 4: (color online) (a) Temperature dependence of the 
block AFM moment at n = 2.0. Black and red curves are for 
next nearest neighbor s±-wave and nearest neighbor d-wave 
pair couplings, respectively, (b): the re-plotted curve of the 
neutron data from Ref. [ill . 



ment of the BAFM order highly depends on the coupling 
constant Vr^ij. For realistic parameters of the interlayer 
tunneling, our calculation shows that the superconduc- 
tivity proximity effect may result in substantial suppres- 
sion of the magnetic moment. This is in contrary to that 
in the cuprate superconductor, where the coupling con- 
stant Vr,ij is very small because of small c-axis hopping 
integral due to the large anisotropy, and because of the 
renormalization of Vr^ij by a factor proportional to hole 
concentrations due to the no double-occupation condi- 
tioni^ In iron-based superconductor, the anisotropy of 
iron-based material is much smaller than in the cuprate, 
which lead to a relative larger t^-- And the moderate 
correlation effect in iron-based superconductor leads to a 
moderate renormalization factors. As a consequence, the 
coupling constant VT,ij in iron chalcogenide superconduc- 
tor should be moderate. 

We remark that we'd be careful in drawing a concrete 
conclusion to compare with the experiments. The ap- 
proximation that only dxz and dyz orbitals are impor- 
tant in the bands close to Fermi energy is good in terms 
of the band structures. But the maximum magnetic 
moment in two-orbital model is only 2/iB, smaller than 
the moment of 3.31/LtB measured in experiments i ^ ^ i The 
other effect is that we only calculated the suppression 
of the BAFM order of the surface layer of the BAFM 
domain. According to TEM experiment^ each BAFM 
domain has around ten layers. And the suppression of 
BAFM order of the layers in the middle of domain may 
be more complicated. In brief, the suppression of the 
BAFM moment is sizable because of the moderate cou- 
pling constant Vr,ij, and our calculation may be viewed 
as a semi-quantitative result. 

We also investigated proximity effect for various pair- 
ing symmetry of the SC phase. It has shown that the SC 
pairing with NN d-wave symmetry resulted a more pro- 
nounced proximity effect in reducing the moment of the 
BAFM order than the NNN s±-wave pairing. The second 
order process induced proximity effect has a temperature 
dependent as the SC pairing, which may be relevant to 
the experimental observations. More remarkable proxim- 
ity effect was found in the BAFM state by comparison 
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with the conventional AFM state, which was the con- 
sequence of the frustrated structure and the associated 
anisotropic exchange interactions. 
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VI. APPENDIX 

In the following, we compare the above proximity effect 
with that in the single band conventional AFM (CAFM) 
system. In order to make the comparison more con- 
vincing, we choose the dispersion Sk = —2t[cos{kx) + 
cos(fcy)] — At' cos{kx) cos{ky) — fi with t = ti and t' — t^, 
which gives rise to the similar energy band width with 
that in the above two-orbital model and is close to the 
case of the cuprates. The AFM order is introduced by the 
AFM exchange interaction JJ2{ij) SrSj between the NN 
sites. At the half filling n = 1, we find that J = 1.6 pro- 
duces the comparable band gap and the electron polariza- 
tion as in the above BAFM state. In Fig. [5l we present 
the magnitude of the magnetic order and the induced 
pairing correlation as a function of the effective tunneling 
Vr^ij. The upper panel shows the results for the NNN s±- 
wave pairing and the lower panel the results for the NN 
d-wave pairing. In the figure, the magnitude of the mag- 
netic order in both cases is renormalized. The proximity 
effect in reducing the AFM order is more pronounced for 
the BAFM state as shown in Figs. [5ja) andlSjc). As for 
the induced pairing correlation, the larger correlation is 
found in the BAFM state for the s±-wave paring and in 
the CAFM state for the d-wave pairing, as displayed in 
Figs.[5jb) and[5jd), respectively. 

We can understand the above results by considering 
the different spin configurations of the BAFM and CAFM 
orders, as shown in Fig.[6l In the BAFM state, when two 
electrons transfer from the BAFM layer to the SC one, 
the energy changes due to the bonds breaking for the 
NN bond coupling are AeH^ = \Ji\ [A1 and A2 bonds 
in Fig. Eta)] and AE^fj^ = 7| Ji|/4 [A3 bond in Fig. Ha)] 
for the electron pairs with ferromagnetic and antiferro- 
magnetic alignments, respectively. On the other hand, 
the energy changes due the bonds breaking for the NNN 
9|Ji|/4 [B3 bond in Fig. Eta)] and 

As 

a result, the proximity effect in reducing the moment of 
the AFM order by the d-wave pairing is more remark- 
able than that by the s± one. However, the energy 
changes due to the bonds breaking in the CAFM state are 



bond are AeJJj^j^ 
AeJ^,^j^ = 5| Ji|/2 [Bl and B2 bonds in Fig. Elja) 
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FIG. 5: (color online) Comparison of the proximity effect 
between the block and conventional AFM states. Left column 
shows the moment of the AFM order, and right column the 
SC pairing correlation as functions of the effective tunneling 
strength Upper panel: for the next nearest neighbor s±- 

wave pairing and lower panel for the nearest neighbor d-wave 
pairing. 
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FIG. 6: (color online) Comparison of the spin structures and 
their respective NN and NNN bonds, (a): block AFM state; 
(b): conventional AFM state. 



AEnn = 7\J\ [C bond in Fig.EIb)] and AEnnn = S\J\ 
[D bond in Fig. Elb)] for the NN and NNN band cou- 
plings. Therefore, the proximity effect in reducing the 
moment of the AFM in the CAFM state is rather weak 
for both the s±- and d-wave pairing couplings. As for 
the induced pairing correlation, the extent of the match 
between the AFM order configuration and the singlet SC 
pairing largely determines the magnitude of the induced 
pairing correlation. For example, the CAFM matches 
well with the NN d-wave pairing, so that one can ex- 
pect a large induced pairing correlation without the se- 
vere decrease of the AFM order as displayed in Figs.[Sl[c) 
and Eld). 
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